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ABSTRACT 



Optical photometry, spectroscopy and XMM-Newton ultraviolet and X-ray observations with full phase coverage are used for an in-depth study 
of WXLMi, a system formerly termed as a low-accretion rate polar We find a constant low mass accretion rate, M ~ 1.5 x IO^'^Mq yr"', a 
peculiar accretion geometry with one spot not being accessible via Roche-lobe overflow, a low temperature of the white dwarf, T^g < 8000 K 
and the secondary likely being Roche-lobe underfilling. All this lends further support to the changed view on WXLMi and related systems 
as detached binaries, i.e. magnetic post-common envelope binaries without significant Roche-lobe overflow in the past. The transfer rate 
determined here is compatible with accretion from a stellar wind. We use cyclotron spectroscopy to determine the accretion geometry and to 
constrain the plasma temperatures. Both, cyclotron spectroscopy and X-ray plasma diagnostics reveal low plasma temperatures below 3 keV on 
both accretion spots. For the low m, high B plasma at the accretion spots in WX LMi, cyclotron cooling is dominating thermal plasma radiation 
in the optical. Optical spectroscopy and X-ray timing reveal atmospheric, chromospheric and coronal activity at the saturation level on the 
dM4.5 secondary star 

Key words, stars: individual: WXLMi- stars: magnetic field - X-rays: stars 



1. Introduction 

Polars are magnetic cataclysmic binaries consisting of a late- 
type main-sequence star and a stro ngly magnetic w hite dwarf 
locked in synchronous rotation (see IWarneii (119951) for a com- 
prehensive survey of cataclysmic variable stars). In normal po- 
lars accretion happens via Roche-lobe overflow and accretion 
streams towards the magnetic poles where the accretion en- 
ergy is released mainly as X-ray thermal radiation, optical cy- 
clotron radiation and a prominent soft X-ray component. The 
soft component makes them prominent sources in the soft X- 
ray sky and they were found numerously in the ROSAT all- 
sky survey. The detection bias is large and the true space den- 
sity highly uncertain. Recently a few systems with very low 
accretion rate (a factor 100 - 1000 below the canonical val- 
ues for Roche-lobe over flow) were uncover ed in o ptical spec 



WXLMi 



Reimers et al 



(HS 1023-H3900) was discovered by 
1999h from the Hamburg Quasar Survey. 



trosc opic surveys (HQS: iHagen et al.lll995L SPSS: lYork et al 



2000). They were identified on the basis of highly peculiar opti- 
cal spectra with broad cyclotron harmonics in emission mimic- 
ing quasar emission lines. They could form an important part 
of the population of close interacting binaries, being either nor- 
mal CVs in extended low states or pre-CVs on their track to- 
wards R oche-lobe overflow. R ecently, pushed forward by the 
work bv lSchmidt et alj (l2005h evidence grows to regard these 
systems as pre-CVs. They might serve as tracers to uncover the 
unbiased sample of magnetic CVs. 



From optical and spectroscopic follow-up observations they 
determined a period of about 167 minutes and two accretion 
spots with a field strength of 60 MG and 68 MG, respec- 
tively. The accretion rate was found to be M < 3 x 10"'^ 
Mo/yr, i.e. orders of magnitudes below the normally observed 
accretion rate in polars. A spectral type of dM4.5 gave a 
good spectral fit for the s econdary and w as used to deduce a 
distance of 140 ± 50 pc. ISchwarz et al. I (Uooi ) used optical 
photometry, performed with the Potsdam 70 cm telescope in 
UBVRI, to determine a stable long-term ephemeris based on 
the timings of optical maxima from the primary accretion spot. 
During the whole monitoring campaign covering half a year 
the system was in a similar low state of accretion with the one 
exception of a flare which was located on the active secondary. 
From the colors at the orbital phase when both accretion spots 
are invisible the spectral type of the secondary was likewise 
determined to be M4.5. Using the relative magnitude in the I 
band at photometric minimum yields a distance of ~ 100 pc. 

WX LMi as well as t he similar objec t s terme d LARPs (Low 



Accretion Rate Polars, ISchwope et alj (120021) ) was not dis 



covered in the ROSAT All Sky Survey despite its relatively 
short distance. WX LMi was in the off-axis field of two X-ray 
ROSAT PSPC pointings (10 and 13 ksec) and discovered at a 
rate of 0.004 cts/sec. The very low number of photons did not 
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Table 1. Log of spectroscopic observations with the Calar Alto 
3.5m telescope 



Date 
(Y/M/D) 


Instrument 


Wavelength 
range [A] 


Resolution 
[A] 


No. of 
spectra* 


1999/3/9 


TWIN 


3800-9900 


6/4° 


47 


1999/3/10 


TWIN 


3800-9900 


6/4° 


56 


1999/3/18 


MOSCA 


3380-8410 


13 


29 



° for the blue/red channel 



* integration time 300 seconds each 

allow to extract a proper X-ray spectrum and to determine the 
origin of the X-rays, being either accretion-induced or from the 
active corona of the secondary. We thus performed new X-ray 
observations with XMM-Newton in order to study magnetic 
accretion in an extreme combination of the parameters, which 
control the process of energy release, i.e. at very low specific 
accretion rate m and at high magnetic field. The X-ray observa- 
tions of WX LMi were accompanied by near-ultraviolet obser- 
vations with the OM onboard XMM-Newton. It was used as a 
bolometer in order to determine the white dwarf photospheric 
temperature. This should help to decide if the LARPs are nor- 
mal polars in an occasional state of very low accretion or pre- 
cataclysmic binaries that never have been accreting. Our newer 
X-ray and UV observations are analysed together with phase- 
resolved optical spectroscopic observations performed already 
back in 1999. Since the overall brightness of WX LMi does not 
change, the two data sets can be combined without problem 
and thus allow a broader picture of the relevant spectral com- 
ponents. 

The paper is organized as follows. In Section |2] the new 
optical and X-ray observations and the reduction steps are de- 
scribed. In Sect. 3 we begin the analysis by an in-depth study 
of the spectral features from the secondary which constrains 
the spectral types, the distance, its activity and the binary pa- 
rameters. In Sect. 3.2 we describe our modeling of the phase- 
resolved cyclotron spectra, which constrains the locations of 
the accretion spots and the temperatures in the accretion plas- 
mas. Sect. 3.3 describes our attempts to understand the OM- 
filter observations and Sect. 3.4 finally gives an analysis of the 
X-ray spectrum. 

2. New optical and X-ray observations 

2.1. Optical photometry - updated ephemeris 

Schwarz et al.l (l200ll) fixed the ephemeris of the system using 
photometric data covering a period of six month. The primary 
maxima of the V band data pinned down the epoch of zero 
phase and the variability of the light curves in VRI was used to 
establish the period. We performed additional photometric ob- 
servation with the Potsdam 70cm telescope in V band on April 
28 2004 and March 21 2005. These observations revealed the 
system to be at the same brightness as during our monitoring 
observations in 1999 and also showed the same variability pat- 
tern. A linear regression to all observed primary V band max- 
ima between 1999 and 2005 yield a period of 0. 1 1592364(43) 
days. Our spectroscopic observations from March 1999 deter- 
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WAVELENGTH (A) 

Fig. 1. Mean orbital spectrum of WX LMi obtained 
March/April, 1999, with the TWIN and MOSCA spectro- 
graphs at the Calar Alto 3.5m telescope. The Ca H-i-K, the 
Balmer emission lines and Na absorption lines are indicated 
by vertical ticks. 

mine the inferior conjunction of the secondary. The epoch of 
the blue-to-red zero crossing of the Nal-lines (see Sect. I3.1.3I I 
was used as the spectroscopic zero point. This yields the up- 
dated ephemeris 

BJD = 2451247.50425(27) -H £ x 0.11592364(43) (1) 

where the number in brackets represent the uncertainty. The 
difference between spectroscopic and photometric phase zero 
is 0.168 phase units. The phase used in this paper refers to this 
spectroscopic ephemeris. 

2.2. Optical spectroscopy 

WXLMi was observed with the Calar Alto 3.5m telescope at 
two occasions in March 1999. The double-beam spectrograph 
TWIN was used for the first run in March 9-11, the multi- 
object spectrograph MOSCA for the second run on March 
18/19. Low-resolution gratings with reciprocal dispersions of 
144 A/mm and 108 A/mm were used for the TWIN obser- 
vations resulting in a spectral coverage of 3800-6900 A and 
6700-9900 A in the blue and red channels at a spectral reso- 
lution of about 6 and 4 A, respectively. A series of 47 spec- 
tra with exposure time of 5 minutes each were obtained in the 
night March 9/10 (BJD 2451247.484406 - 2451247.671117), 
another series of 56 spectra with the same exposure time in the 
following night (BJD 245 1 248 . 375258 - 245 1 248 .596263). All 
binary phases were thus covered by spectroscopic observations. 

Spectrophotometric standards were observed in the same 
nights and used for calibration of the spectral response. 
However, the weather conditions were rather poor during the 
TWIN observations with transparency variations by 100%. 
Hence, we could not use the standard star observations to 
put the spectra of WX LMi on an absolute scale, we could 
just determine the spectral slope. The nearby star GSC2- 
N23233202253 (a (2000): 10:26:27.5, 6 (2000): -h38:45:03, 
distance to WXLMi: 86") was used as comparison star and 
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also placed on the spectrograph's slit. We tried to use the vari- 
able signal of this star for a photometric calibration of the 
spectra of our target. This attempt turned out not to be suc- 
cessful since the photometric variations showed color-terms, 
likely to be different between the target and the comparison 
star (likely/possible reasons for that being a slight misaligne- 
ment of the spectrograph slit, observation not being performed 
at the parallactic angle, the different spectral slopes of the target 
and the comparison). 

Instead we used our photometric observations at the 70cm 
telescope to calibrate our spectra (assuming no overall change 
of brightness and color). This approach was also not found to 
be free of problems since the wavelength ranges of the spectro- 
graph channels were too narrow to cover more than one (red) 
or two (blue) filters completely which made a correction of the 
color terms rather uncertain. 

He/Ar arc lamp spectra were observed for calibrating the 
dispersion. Flexure of the spectrograph resulted in considerable 
shifts of line positions on the CCD of up to 0.6 pixels. These 
were measured in each spectrum individually by correlating 
the night-sky spectrum with a template. These time-dependent 
shifts were then approximated by a low-order polynomial and 
the dispersion relation was shifted per spectrum according to 
those fits. The late-type stars G1213, G1268, and G1273 were 
observed as possible template stars for the main-sequence sec- 
ondary in WX LMi through wide and narrow slits in the night 
of March 10. 

We followed a similar observational strategy for the 
MOSCA observations by placing the same comparison star on 
the spectrograph's slit and obtained 29 exposures of 5 minutes 
each in the night March 18/19 between UT 22:45 and 1:41 
(BJD 451256.454120-2451256.576432), i.e. we covered a bit 
less than one orbital cycle of the 167 min binary. The B500 
grism was used as a disperser resulting in a spectral coverage 
3380 - 8410 A with about 13 A (FWHM) resolution. Similar to 
the TWIN observations, small remaining shifts of the positions 
of the night sky lines in each spectrum were used to improve 
the dispersion relation which was originally based on HgNeAr 
arc lamp spectra. 

Also the MOSCA observations suffered from large trans- 
parency variations by almost 100%. This time these variations 
could rather successfully corrected for with the help of the 
comparison star observed through the same slit. A further com- 
plication arose from the fact that no spectrophotometric stan- 
dard stars could be observed during that night. A photometric 
calibration of the MOSCA spectra was achieved by assuming 
that the orbital mean spectrum of WX LMi was unchanged be- 
tween the TWIN and the MOSCA observations, an assumption 



Table 2. Parameters for the secondary in WX LMi 



which is justified by our long-term photometry ( ISchwarz et al 



20011) . Division of the uncalibrated mean MOSCA spectrum 
by the calibrated mean TWIN spectrum yielded an instrumen- 
tal response curve for the MOSCA observations. 

2.3. XMM-Newton X-ray and ultraviolet observations 

We performed ultraviolet and X-ray observations of WX LMi 
with XMM-Newton in April 2004 for about 30 ksec, cov- 



Sp 


log{R/R^r 






5'' 


' TiO 


d{SK) 


rf'(TiO) 


4 


-0.613 


0.220 


7.53 


4.37 


3.08 


103 


130 


4.5 


-0.676 


0.179 


7.95 


4.46 


2.62 


85 


104 


5 


-0.740 


0.147 


8.37 


4.56 


2.15 


70 


82 



" according to Beuermann et al. ( 1999) 
* according to Delfosse et al. ( 2OO0 
^ according to Delfosse et al. (2000 
'' according t olBeuermann et al.l jl999.. 
' according to lBeuermann et aL Ul999l) in units of 10 erg cm s A 

using an observed flux deficit F7165 = 5.5 x 10"'* erg cm" 
(uncertainty of 30%) 



ering 3 orbital cycles. Calibrated photon event tables for all 
detectors were computed with SAS version 6.5. Since the 
spectrum turned out to be very soft, the SAS task epreject 
was used to extend the usable energy range down to 0.12 
keV. Barycentric correction was applied using barycen. The 
spectra were extracted from the event tables with SAS tasks 
evselect and especget. The light curve was obtained run- 
ning the edetect_chain pipeline on time selected event tables 
containing all events ocurred in a given orbital phase interval. 

The EPIC PN was operated in full frame imaging mode, 
yielding ~ 800 counts from the source. Both EPIC MOS cam- 
eras were operated in partial window mode, yielding both ~ 
200 counts. The observing time with the OM was spent in 
nearly equal shares with the U and UVWl filters centered on 
3440 A and 2910 A, respectively. The detector was used in fast 
mode for the full observation. The mean countrates in the U 
and UVWl filters were 2.12 s"'and 0.75 s"', respectively. 

3. Analysis 

Fig. [U shows the mean orbital spectrum of WXLMi ob- 
tained with the TWIN (longward 6800 A) and MOSCA 
spectrographs . It sh ows the same features as described by 
Reimers et al. d 1999b . .e. an M dwarf dominating the red spec- 
tral range, a white dwarf which is responsible for the spec- 
tral upturn in the blue range and several pronounced cyclotron 
lines, the most prominent one at 5950 A. The higher spectral 
resolution and the partly better signal provided by the observa- 
tions presented here allows individual spectral features to be re- 
solved and studied. The most prominent ones are the H-Balmer 
emission lines and the Nal absorption lines at 8 183/94 A. We 
therefore begin our analysis with an optical spectroscopic study 
focusing on the secondary star and the optical cyclotron spec- 
trum. 



3.1. The secondary in WXLMi - optical spectroscopy 
3.1 .1 . Spectral type of the secondary 

In order to determine the spectral type of the secondary we 
make use of the TWIN spectra obtained in the red channel. 
The spectral type is not determined straightforwardly due to 
the wavelength- and phase-dependent background component 
from the white dwarf (photospheric and beamed cyclotron ra- 
diation from two accretion spots). The remaining calibration 



4 



J. Vogel, A.D. Schwope and B.T.Gansicke: An in-depth study of the pre-polar candidate WXLMi 



uncertainties, as described above, gave rise to some additional 
complication. 

We firstly shifted all spectra to radial velocity zero using the 
sine fit to the Na radial velocity curve (see below). We then se- 
lected only those spectra which were not too heavily affected 
by atmospheric absorption, i.e. those with the best signal-to- 
noise ratio, and calculated an average spectrum. 
We adapted a non-magnetic model atmosphere of the white 
dwarf to the average spectrum in a wavelength region free of 
any cyclotron emission and subtracted the same scaled model 
spectrum from all selected observed spectra. This left us with 
phase-dependent spectra containing only the secondary and the 
cyclotron lines. 

In order to constrain the spectral type of the secondary we de- 
cided to use the narrow-band spectral indices Ti05 and VO-a 
jCruz & Reid (2002), Kirkpatrick et al. (1999)). According to 
those papers these indices cor relate best with the spectral type. 
We also tried the PC3 index (" Martinet al]ll999 ') although we 
did regard the results less reliable since the two continuum side 
bands have a relatively large separation and the index is there- 
fore much more sensitive to variations of the underlying cy- 
clotron radiation. 

We determined the indices for all spectra, eliminated those in- 
dices where the spectral features were found superposed on 
a cyclotron line and finally calculated average indices. The 
Ti05 index thus determined was 0.30 + 0.05, the VO-a in- 
dex 1.99 + 0.03, and the PC3 index to 1.21 + 0.05. The given 
errors reflect the statistical uncertainty and the uncertainty 
which originates from the subtraction of the white dwarf model 
spectrum. The measured indices resulted in spectral types of 
M5.0+0.5, M4.5+0.7, and M4.5+1. We finally adopted a spec- 
tral type of M4.5 with an uncertainty of half a subclass. This 
new determination of the spectral type based on spectral data 
with su fficient resoluti o n is in good agree ment with th e esti- 
mates of iReimers et al ] (Il999h and lSchwarz et al.1 (1200 ih based 
on spectroscopy with much lower resolution and broad-band 
optical pho tometry, respectively. 

Beuermann et al. I d 19981) have shown that the color I - K can 
be used as tracer of the spectral type. A parameterization of 
the Sp-(/ - K) relation was given as a third-order polyno- 
mial. Due to the high field strenght in WXLMi the /iT-Band 
is likely uncontaminated by cyclotron radiation (see Fig.[T2lfor 
the wavelength of the cyclotron fundamental for both poles), 
thus allowing to dete rmine a spectral type of M4 for / - 14.83 
dSchwarz et al .11200 i f) and K = 1 2.49 (2MASS d atabase) in ac- 
cordance with the above value. iBeuermann et al. (1998) have 
also shown that most CV secondaries with Po^b between 2 and 
5 hours are cooler than expected for Roche-lobe filling ZAMS 
stars. Nuclear evolution prior to mass transfer and lack of ther- 
mal equilibrium due to mass loss were mentioned as possible 
causes of this discrepancy. We will argue below that both op- 
tions likely do not apply here suggesting the presence of a 
Roche-lobe underfilling main-sequence secondary in WX LMi. 




VELOCITY (km/s) 



-2000 2000 

VELOCITY (km/s) 



Fig. 2. TWIN spectroscopy of WXLMi. Shown are cutouts 
of the phase-averaged spectra arranged as an apparent trailed 
spectrogram centered on the Ha emission line and the Nal ab- 
sorption lines. Phase zero refers to inferior conjuntion of the 
secondary, Doppler shifts were transformed to velocity. 



3.1 .2. The distance to WX LMi 

The composite nature of the emission from cataclysmic vari- 
ables, particular in the optical, makes it difficult to perform a 
straightforward distance measurement. We do not make use of 
Roche lobe geometry, since the secondary is possibly some- 
what underfilling. The simple use of the distance modulus also 
doesn't seem to be appropriate, since the secondary is clearly 
deformed as evidenced from ellipsoidal modulations (Fig. [H 
bottom panel). We assume that the secondary in WXLMi is 
a ZAMS star with solar metallicity and make use of the sur- 
face brightne ss - spectral type and Fjio - spectr al type relations 
compiled by IBeuermann & Weichhold (Il999l). We use radii 
from the compilations of Beuermann et alj (Il999 *). Absolute 
magnitude s and masses for the p ossible spectral types are de- 
rived from lOelfosse et al.l (120001) . The results are summarized 
in Table |2l 

The distance via the surface brightness S k assumes that the 
2MASS /T-band magnitude is solely due to the secondary star 
The measured flux deficit in the TiO feature was corrected for 
the contribution of the white dwarf at this wavelength. For the 
most likely spectral type, M4.5, our distance estimates give 
85 and 100 pc, respectively, with a relatively large spread of 
about 20 pc. We regard the derived distance via the S k relation 
somewhat more rehable, since the photometric data are not af- 
fected by the calibration problems of our spectroscopy. If the 
secondary would behave like most of the CV secondaries at 
the given Poib, the radii would be larger than given in Tab. |2] 
and the distance would increase. We thus use in the following 
a distance of lOOpc to WX LMi. 



3.1.3. Emission and absorption line radial velocity 
variations 

Our final reduced spectra were phase-averaged and arranged 
as an apparent trailed spectrogram using our improved orbital 
period and the phase of inferior conjunction of the secondary as 
determined below. A cutout of the spectra centered on the Wa 
emission lines and the Nal absorption lines is shown in Fig. |2] 
Wavelength was transformed to radial velocity in the diagrams. 
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Fig. 3. Radial velocity of the secondary as a function of the 
mass ratio and the orbital inclination (between 50° and 90° in 
steps of 10° ). The horizontal dashed lines mark the observed 
radial velocity with Icr uncertainties. 

Table 3. Radial velocity measurements of WXLMi in 1999. 
Phase zero is defined by the spectroscopic ephemeris at BJD 
(0 = 0) = 2451247.50425(27) + Ex 0.11592364(43) 



Line 



Date 



K (km/s) 



00 



Nal (TWIN) 
Ha (TWIN) 
Ha (MOSCA) 



March 9-1 1 
March 9-1 1 
March 18/19 



286: 
330: 

268 : 



0.0! 
0.008 + 0.002 
0.005 ± 0.004 



Nal-line flux at phase 0.5 with respect to phase 0.0. However, 
the flux measurements of the Nal-Unes are rather uncertain due 
to a too low signal-to-noise ratio of the individual spectra and 
we didn't analyse this line property in more detail. 

The tight phase relation between the Ha and the Nal lines 
suggests a common origin on the secondary star This is sup- 
ported by the absence of any further emission line component 
in Ha which could be attributed to an accretion stream or disk. 
The Ha emission line is assigned completely to the photo- 
spheric activity of the secondary star 

We regard the radial velocity amplitude K of the Nal lines 
as proper measurements of the projected orbital velocity of the 
secondary. Irradiation effects, which distort the radial velocity 
curves in high-accretion rate polars, play an insignificant role 
in WXLMi. The observed velocity amplitudes of Ha in the 
TWIN and MOSCA spectra are inconsistent with each other 
and with Nal. Hence, the centers of activity on the secondary 
were located away from the center of mass and were time vari- 
able on a scale of a few days. 

The observed radial velocity amplitude gives some con- 
straints on the mass ratio and the mass of the white dwarf 
(see Fig. O. The absence of an eclipse constrains the orbital 
inclination to ; < 72° assuming a Roche-lobe filling secondary. 
For a Roche-lobe underfilling star the inclination could be even 
higher With the ZAMS assumption for the secondary and the 
corresponding masses as listed in Tab.|2]the minimum mass of 
the white dwarf is 0.36 Mq. 



The lines from both species are clearly detected through the 
complete orbital cycle. 

The Nal-doublet is clearly resolved in the TWIN spectra. 
Radial velocity measurements were performed by fitting dou- 
ble Gaussians with fixed wavelength separation to the observed 
line profiles. Those attempts were not always successful due to 
the low signal to noise ratio in several of the spectra. For the 
determination of the radial velocity curve we neglected spectra 
below a certain S/N level and were thus left with a total of 65 
spectra. 

Despite the lack of sufficient resolution the Ha lines of 
WXLMi seem to consist of just one line. Contrary to polars 
in their high accretion states, which show even at our spec- 
tral resolution more complex line profiles. Ha line positions 
in the individual spectra were thus determined by fitting sin- 
gle Gaussians. The Ha lines were detected also in all MOSCA 
spectra and line positions were determined by fitting Gaussians 
in 26 spectra. 

The radial velocity curves of Nal and of Ha at both oc- 
casions could be successfully fitted with sine curves. A sum- 
mary of the radial velocity measurments is given in Table [3] 
The epoch of the blue-to-red zero crossing of the Nal-lines at 
BJD = 2451247.50425(27) determines inferior conjunction of 
the secondary and spectroscopic phase zero throughout this pa- 
per. The number in parentheses gives the uncertainty in the last 
digits of the ephemeris zero point. 

We could not discover any significant variation of the width 
of the Nal-lines. There is a tendency of a 20% increase of the 



3.1 .4. Emission line fluxes as activity indicators 

The Ha and Ca H&K emission lines are regarded as originat- 
ing from the active photosphere and chromosphere of the sec- 
ondary. We arranged our observed spectra like a trailed spec- 
trogram by phase-averaging the orginal spectra into 30 phase 
bins. We then shifted the spectra to radial velocity zero and sub- 
tracted the scaled spectrum of the single, non-active M4.5 star 
Gl 268. The remaining flux in the spectral region around Ha 
consisted of the smoothly varying continuum from the white 
dwarf and the Une flux without Doppler shift and thus allowed a 
straightforward measurement of the line flux. The Ha emission 
line flux was related to the continuum of the scaled Gliese star 
at this wavelength in order to determine the equivalent width of 
the line. The results of this procedure are displayed in Fig.|4] 

The continuum emission light curve shows a double- 
humped structure with minima at phase and 0.5, respectively, 
very reminiscent of ellipsoidal light variation. The /?-band light 
curve of the scaled template spectrum shows modulations of 
0.2 - 0.3 mag. Using nightfallQ we made some attempts to 
model this pattern. At a temperature of 3300 K for Sp M4 one 
needs a Roche-lobe filling secondary and a large inclination an- 
gle of order 70° in order to get such a big amplitude. However, 
the photometric accuracy achieved by us does not allow to draw 
firm conclusions on the filling factor of the secondary, on the 
temperature or the inclination. E.g. the maxima of the contin- 



' http : // www ■ hs ■ uni -Hamburg . de/DE/Ins/Per /Wi chmarm/ 1 
Nightfall.html 
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Fig. 4. Phase-dependent Ha emission line and continuum flux 
as well as the equivalent width as measured in the TWIN spec- 
tra. 



uum light curve are differing by ~0. 15 mag, probably a leftover 
from our calibration uncertainties. 

The line flux itself does not show a clear phase-dependent 
variability pattern, whereas the eqivalent width again shows a 
pronounced doubled-humped structure with maxima at phase 
and 0.5, respectively. The EW variied between 4 and 8 A. The 
subtraction of the template M-star spectrum had some uncer- 
tainty, using a possible different scaling raised EW up to 10 A. 
These numbers were clearly smaller during the MOSCA obser- 
vations with a variation through the binary orbit between 1 and 
4 A. 



WaUcowicz et al.l ( 120041) investigated the strength of activity 



via Hot measurements for a large number of SDSS stars. They 
parametrise the strength with ax factor, x EWh„ = Ln^/Lboi- 
For spectral types M4 - M4.5 the;^' factor is log;^' = -4.2 . . . - 
4.4 and log(LH„/iboi) becomes -3.7 ... - 3.2 for the TWIN 
observations. 



Recently. I West et al.1 (12004ft and Silvestri et al.l (120051) stud 
ied large number of either single (IWest et al .1 12004 or com- 



mon p roper motion WD/MS binaries (CPMBs, ISilvestri et al 
(12005b') drawn from the SDSS. The large majority of 90% of 
the lSilvestri et all ( l2005b CPMBs have an EW close to (their 
Fig. 2). The mean Ha luminosity of the West et al. sample 
at spectral type M4-M5 is -3.57 ... - 3.67 (their Fig. 4 and 
Tab. 1). These figures clearly put the secondary in WXLMi 
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Fig. 5. Contour plot of the phase-dependent cyclotron spectra 
from the two accretion regions on the white dwarf in WX LMi. 
Frequency is given in units of 10'^ Hz. Numerals at the top indi- 
cate frequencies for the different cyclotron harmonics. Roman 
numerals indicate primary and secondary spot, Arabic numer- 
als indicate harmonic numbers. 



among the more active stars in its category. The comparison be- 
tween our MOSCA and TWIN observations showed the Ha ac- 
tivity remarkably variable on a timescale of a month or shorter. 

The Ca II H&K lines appear as distinct features only in 
the radial-velocity corrected average spectrum of WXLMi 
(Fig. [TJ. They can only hardly be recognised in the original 
spectra where they are immersed in the noise. A continuum un- 
der the emission lines was determined interactively on a graph- 
ics screen and the excess flux regarded as Ca ii H&K emis- 
sion line flux (with a negligible contamination from hydro- 
gen He). This gave an integrated flux of Fqh - (1.5 - 2.0) x 
10"'^ erg cm"^ s"' and a lumino sity of log Lca(erg/s ) ~ 27.3 for 
an assumed distance of lOOpc. iPiters et al.l (Il997h have stud- 
ied the relation between coronal and chromospheric emission 
from cool stars in near-simultaneous ROS AT all-sky survey and 
Mount Wilson data and found a good correlation between these 
two activity indicators. Would the secondary follow this re- 
lation (their Fig. 4), we would expect an X-ray luminosity of 
log Lx{erg/s) ~ 33, far higher than observed (see below). 

Vice versa, the observed X-ray lumonisity of WX LMi im- 
plies a Ca ii luminosity 4 orders of magnitude less than ob- 
served, would the star follow the Piters et al. relation. These 
estimates and comparisons suggest, that either the Lx - Lca 
relation between coronal and chromospheric luminosities does 
not apply to the case of the fastly rotating M star in WX LMi 
or a very high variability of the activity parameters of WX LMi 
or a combination of both. The implied X-ray luminosity based 
on the observed Ca ii luminosity would be higher than the 
typical accretion-induced luminosity in a high accretion polar. 
Probably we have encountered WXLMi in a state of unusual 
high chromospheric activity during the MOSCA observations. 



3.2. Optical cyclotron radiation 

The cyclotron spectra from the two accretion regions on the 
white dwarf in WX LMi were extracted from the original data 
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Fig. 6. Light curves of the 3"* and 4* (thick Hne) cyclotron har- 
monics from the primary (bottom panel) and the 2"'' (thick line) 
and 3"^ harmonics from the secondary (top panel) accretion 
spot in WXLMi. The light curve for the 2"'^ harmonic of the 
second spot is smoothed to reduce the noise level. 



in a similar manner as described by iReimers et al. I (ll999h . We 
firstly corrected the spectra to radial velocity zero, arranged the 
spectra as a trailed spectrogram by phase-averaging the spectra 
into 30 phase bins of same length (this left 4 bins unpopulated 
in the MOSCA data), and subtracted suitably scaled spectra 
of our template Gl 268 (M4.5). This gave us spectra contain- 
ing mainly photospheric radiation from the white dwarf and 
cyclotron radiation from the accretion spots (see Fig. [10] for a 
graph of the spectral decomposition). In the following we are 
using the terms 'accretion region' and 'pole' synonymously, 
just for brevity since we know about the differences. 

Assuming a flat light curve from the rotating white dwarf, 
a mean white dwarf spectrum was composed from spectral re- 
gions free of cyclotron radiation and subtracted from all phase- 
binned spectra. Those diff'erence spectra are regarded as of 
pure cyclotron origin. They were arranged as a trailed spectro- 
gram and are shown on a logarithmic intensity scale in Fig.|5] 
Actually, the described procedure was performed for the red 
and blue TWIN spectra and the MOSCA separately and the 
trailed spectrogram shown in Fig. |5] is a composite of those. 
The MOSCA data are used for the blue end of the spectra, 
V > 6.3 X 10'"* Hz, since their signal-to-noise ratio was supe- 
rior to the TWIN spectra in that wavelength region. The feature 
at V = 4.3 X 10'"^ Hz (0 ~ 0.65) is a leftover from the M star 
subtraction. 

Fig. [5] may be co mpared with the corresponding Fig. 3 of 
Reimers et al.l(ll999h . As a novel feature we discover the 4* cy- 
clotron harmonic from the primary pole at v = 6.72 x 10'** Hz 
in the MOSCA data. All other cyclotron harmonics appear 
ve ry similar in s t rength , shape and position to those shown 
by iReimers et al.l (Il999h . Light curves for individual harmon- 
ics were computed by averaging the trailed spectrogram over 
certain wavelength ranges. The light curves for the 3"^ and 4* 
harmonics from the brighter primary and the 2"*^ and 3"* har- 



monics from the fainter secondary pole are shown in Fig. |6l 
Since the continuum subtraction for the lowest frequencies was 
not completely satisfying, the light curve for the 2"'' harmonic 
of the primary pole is ignored for the detailed analysis. The 
light curves from both poles are double-humped for the higher 
harmonics, an effect of cyclotron beaming in an optically thin 
plasma. The light curves of the lower harmonics are single- 
humped with maximum brightness centered between the two 
beamed humps from the higher harmonics, indicating consid- 
erable optical depth in those harmonics. From the centers of 
the bright phases, 0ci - 0.84 and (p^2 - 0.27, the azimuth of 
both accretion spots were determined to i^i ^ 55° . . . 60° and 
1^2 ^ -90° ... - 100°. 

We modeled the phase-dependent spectra assum- 
ing h omogeneous, isothermal conditions in the plasma 



dChanm ugam & Dulk 1981). A number of parameters deter- 
mine the model spectra: the field strength B, the temperature 
kT, the plasma density coded with the density parameter 
log A, the orbital inclination ;, the co-latitude of the field in 
the accretion spot fij, and the azimuth i/// of the field in the 
spot. The visibility of an emission region as a function of 
phase is further determined by the latitude /3s and azimuth i/Zs 
of the accretion spot and the vertical extent of the emission 
region. For simplicity, we assumed the field in the spot being 
normal to the surface of the star, i.e. we set jS/ - Ps and 

- ips. We further neglected any vertical extent but allowed 
the point-like emission region raised to some height h above 
the star's surface. The orbital inclination was fixed at 70°. The 
difference in azimuth between the two regions was determined 
from the bright phase centers to 150° . . . 160°. 

With these simplifications, we cannot expect to reach a 
perfect fit, particularly not to the light curves of individ- 
ual harmonics which are shaped by the geometry of the ac- 
cretion region, the viewing geometry, optical depth and cy- 
clotron beaming. However, the main features of both cyclotron 
line systems as a function of phase could be reproduced (see 
Fig. |7]i by our model with the following set of parameters: 
Bi,2 = 61.4/69.6MG, logAi^z = 2/2, kTi^2 < 3keV, /!i,2 = 
0.05/0.1 Rwd, j0i,2 = 145°/135° with typical uncertainties of 
AB = 0.5 MG, AlogA = 0.5, AkT = 1 keV, A/i = 0.05 Rwd, 
and A/? = 15°. Some of the parameters are coupled to others, 
e.g. if one raises the temperature one has to lower the density 
parameter in order to match the observed turnover from opti- 
cally thick to thin radiation. The values given here a re slightly 
different from those given in ISchwope et al.l (l2002h based on 
an analysis of the Reimers et al. data. The differences reflect 
mainly the modeling uncertainties and not a difference in the 
data. These were found to be consistent with each other The 
new data provide as addiditonal constraint the beamed fourth 
harmonic from the second pole. 

In our modeling, the co-latitude fi was chosen according to 
the observed phase separation, of the two beamed humps 
in either the 4* (prime pole) or the 3"^ harmonic (secondary 
pole). The value of [5 is related to the orbital inclination / and the 
length of the self-eclipse of an accretion spot, Ac/i, = 1 - A^^, 
via cos(;7rA(/>j) = - cot(/) cot(/?). With the values of p and / fixed 
we had to assume a certain height of the region in order to re- 
flect the length of the bright phase which lasts longer than the 
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Fig. 7. Comparison of the phase-dependent cyclotron model 
(background image) with the observed data (contour lines, the 
same as in Fig.|5]l 



phase separation of the beams. The more extended visib ility 
of the individual spots was used by ISchwarz et al. ( 200l l) for 
their estimate of the accretion geometry. They determined the 
colatitude fi of both spots as 90° - 110° from light curve mod- 
eling for a point source on the surface on the WD. Using the 
beaming proporties of the 3"* harmonic from the second pole 
the range of possible values for the colatitude of the second 
pole was extended up to 130°, which is in agreement with our 
results. Since we are using solely the beaming properties of the 
cyclotron radiation for both spots our results seem to give more 
reliable constraints than the assumption of simple geometrical 
foreshortening. 

These results are unusual as far as the spot locations and 
the plasma temperatures and densities are concerned. The two 
spots are both located in the 'souther n' hemisphere of the white 
dwarf, a result akeady mentioned by lSchwarz et al. ( 200 lb . The 
opening angle between the two spots is of order 160°, i.e. the 
two spots could be located at the two footpoints of the same 
dipolar fieldline. If so, the southern latitude of both spots hints 
to an ofF-centered field configuration. The slight difference be- 
tween the field strengths in the two spots seems to underline 
this, although the strong radial dependence of the field strength 
does not allow to draw firm conclusions about the field at the 
proper footpoints at zero height. One pole lies in the same sec- 
tor where most pole s in high-accretion rate polars are found 
(see ICropperl (119881) Fig. 2), i.e. the sector leading the sec- 
ondary in phase. The other pole is away by 90° from the line 
connecting both stars and away from the direction of a hypo- 
thetical accretion stream. 

The temperatures cannot be determined with very high ac- 
curacy. They are constrained by the width of the observed 
highest harmonic. At temperatures higher than about 3keV 
the predicted line width becomes larger than observed even 
for point-like emission regions. Since in nature there will be 
a spread in kT and B, we regard -3 keV as strict upper limit. 
Thus the temperatures found by us are among the lowest found 
in cyclotron spectra of polars (if WX LMi may be termed as 
such). They are more than an order of magnitude below the 
shock jump temperature for an assumed 0.6 white dwarf. 
The low temperatures together with the low density parame- 
ters A imply that both cyclotron emission regions are domi- 
nated by cyclotron cooling. It seems unlikely, that an accretion 



shock exists in WXLMi, the accretion spots are more likely 
heated by particle bombardm ent. The bombardment models by 
Fischer & BeuermannI (120011) relate the accretion rate per unit 
area, m, to the field strength and the maximum electron temper- 
ature. According to their Fig. 5, one needs for an 0.6 Mq white 
dwarf with rmax/Tshock ~ 0.1 and B - 65 MG a specific mass 
accretion m ~ 0.1 gcm"^ s^'in order to be deeply in the bom- 
bardement regime. 

The mean integrated cyclotron flux calculated from the 
observable parts of the cyclotron spectra are ^ 1 x 
10"'^ erg cm"- s"' and F2 - 0.4 x 10^'^ergcm"^ s"^ re- 
spectively. The correcting factors k for the unobserved parts 
of the cyclotron spectrum and for the beaming pattern are 
rather uncertain, we assume /c = 1-2. We further assume, 
that the cyclotron luminosity can be calculcated as L^yc - 
KnFcfi and get for the sum of both poles Lcyc - (0.4 - 1) x 
10^° ergs~'(c//100pc)^. The implied mass accretion rate of M ^ 
1 .5 X 10"'"^ M0 yr"' is extremely low and is of the same order 
as the wind mass loss rate of the secondary star. 

Given the total and the specific mass accretion rates, one 
may ask for the minimum area over which accretion may hap- 
pen in order to be in the bombardement regime. Equating 
m - M / faccA^Rlj^ and using rii - O.lgcm^^s"', M ^ 
1.5 X lO^'^Moyr-', and 7?wd = 8 x 10'^ cm we obtain for two 
equally sized circular accretion spots a radius of just 40 km or 

/acci— 1.2x10-5. 

3.3. XMM-Newton X-ray and ultraviolet observations 

The main results of the observations with XMM-Newton are 
summarized in Figs.lSlandfTTI Fig.[8]shows the PN X-ray light 
curve as well as the light curves for both OM filters. Along 
with the new satellite data the summarized flux of the 3'^'^ and 
4* harmonics from the primary accretion spot and the 2"'' and 
3'^'' harmonics from the secondary accretion spot is shown. The 
data are binned into 20 phase bins. Around phase 0.6 when 
both accretion spots are self-eclipsed, the cyclotron flux drops 
to zero. The OM and PN light curves exhibit a clear variability 
correlated with the cyclotron flux, and thus the visibility of the 
accretion spots on the white dwarf, which is analysed in more 
detail in the following section. 

3.3.1 . UV - the white dwarf atmosphere 

The UV-observations with the OM were designed to estimate 
the temperature of the white dwarf and of the heated photo- 
sphere below the accretion spots. The basic assumption is that 
the UV light is completely dominated by photospheric radi- 
ation from the white dwarf's surface. Any variability is thus 
assigned to temperature inhomogeneities of the atmosphere. 
We spent each half the available observation time for the U 
(3440A) and UVWl (2910A) filters. We detected a clear, 
phase-dependent variability of the OM count rate in both filters 
correlated with the visibility of the accretion spots (see Fig.[8]i. 
The ratio of the count rates between the two filters, however, 
was constant within the errors indicating only marginal T- 
variations. 
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Fig. 8. Phase folded light curves obtained with EPIC-PN and 
OM together with the total cyclotron flux from Fig. |6]in units 
of 10""' erg cnT^ s"'. 



White dwarf model spectra for pure H, \ogg - 8, non- 
magnetic atmospheres in the temperature range 8000 K - 
1 00000 K jGansickeet al.lll995h were folded through the re- 
sponse curves of the OM with the different filters thus predict- 
ing a count rate in the given filter The count rate ratio between 
the two filters as a function of temperature is shown in Fig.|9] 
The ratio is a sensitive function of T^ff below 30000 K. In order 
to test the reliability of our 'bolometer' we retrieved calibration 
observations of the well-studied white dwarfs BPM 16274 and 
GD 153 from the XMM-Newton Science Archive. Their count 
rate ratios are also shown in the Figure. 

The resulting count rate ratio predicts a temperature of 
39000!^™° K for GD 153 in agreement with values found 



in references listed in the white dwarf data base|j|, 37900 - 
40100K. The derived temperature of 17340!if/*K for BPM 
16274 is somewhat lower than the values of 18400 - 18700K 



15000 



T [K] 



Fig. 9. Theoretical count rate ratio U/UVWl using non- 
magnetic model spectra. Overplotted are measured count 
rate r atios for BPM 16274 with a temperature of 18700 K 
( Bra gagliaetalJl995 ) and GD 1 53 with a temperature of 39600 
K ( Ve nnes et al.ll997h . Error bars have similar sizes as the sym- 
bols. 



found in the literature. We thus estimate the accuracy of the 
filter ratio method of order 10%. 

The observed count rate ratio around phase (p - 0.6, when 
both spots are hidden from the observer, is 2.88 + 0.26. This 
high ratio cannot be reproduced by our model atmospheres. 
The observed ratio indicates a softer spectrum, i.e. cooler atmo- 
sphere than the lower limit of 8000 K of our models. Synthetic 
magnetic white dwarf model spectra with dipolar field struc- 
ture for a pole field strength of 60 MG - kindly provided by 
S. Jordan - reduced this discrepancy, but also led to no satisfy- 
ing result. It need to be said, however, that the model spectrum 
of a magnetised atmosphere used the temperature structure of 
a non-magnetic atmosphere, hence we cannot quantify to what 
extent the magnetic model matches the observations better than 
the non-magnetic. 

It is worth noting at this stage, that we could not identify 
any depression (a line or a trough) in the optical continuum 
of the observed white dwarf spectrum, which could be even 
only tentatively identified as a Zeeman signal. We guess, that 
the low Teff and/or a complex field structure are responsible for 
this observation. 

We conclude that Teff of the white dwarf is almost certainly 
below 8000 K. Even without using the count rate ratio as an 
indicator for the temperature the measured count rate in the 
UVWl filter around phase - 0.6 gives an upper limit for 
the flux and thus the temperature of the white dwarf. Therefore 
one has to assume that the whole white dwarf is seen, with- 
out contribution from the heated pole caps. If we fix the dis- 
tance at 100 pc as derived in section U.1.2l a nd use a mean po- 
lar w hite dwarf mass of 0.6 Mq ( Wickramasinghe & Ferrario 
20001) . the lower limit of 8000 K of our nonmagnetic model 



http://procyon.lpl.arizona.edu/WD/ 



atmospheres gives a flux which is about 12% higher than the 
measured flu x in the UVWl fil ter. We use the mass-radius- 
relation from Nauenber^ ( 1972 ). which is sufficient for very 
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Fig. 10. Observed and deconvolved optical spectrum and ultra- 
violet photometry of WX LMi. Shown are the observed mean 
spectrum (black top curve), the scaled spectrum of G1268 (red), 
the M-star subtracted spectrum (cyclotron radiation plus white 
dwarf, blue). The rhombs denote the OM photometry and the 
black lines at bottom white dwarf model spectra for a distance 
of lOOpc. The black dotted line is the non-magnetic white 
dwarf model with M^d = 0.6 M©, the black solid line denotes a 
magnetic white dwarf model with Mwd - 0.47 Mq, B = 60 MG, 
d = lOOpc. 



cool white dwarfs JPanei et alllioool) . Any additional contribu- 
tion from the accretion spots decreases the flux contributed by 
the white dwarf and thus its temperature. A significantly higher 
temperature would only be possible for a very massive white 
dwarf or if it were much further away than the derived distance 
of 100 pc. 

The different spectral components in the optical and near 
ultraviolet are illustrated in Fig. [TO] which shows the mean 
optical spectrum, the scaled spectrum of the M4.5 template 
(G1268), the blue TWIN spectrum after subtraction of the 
M dwarf (white dwarf plus cyclotron radiation), the two UV 
measurements and suitably scaled magnetic and non-magnetic 
wh ite dwarf mode l specfr a for Jeff = 8000 K. 

Reimers et al. ( 19991) determined the temperature of the 
white dwarf in WXLMi to be 13000 + lOOOK. Our new de- 
termination differs significantly. Reimers et al. substracted a 
dM4.5e (G3-33) spectrum from a mean spectrum when both 
spots are self-eclipsed. The spectral energy distribution of the 
remaining spectrum was then compared with synthetic spec- 
tra for white dwarfs with a magnetic field of 60 MG. The 
temperature dependence of the continuum slope in the wave- 
length range between 4000 to 6000 A used by them is not very 
pronounced. Due to the much smaller contribution of the sec- 
ondary to the total flux in the UV, the use of the OM data seems 
more reliable to us for a determination of the white dwarf's 
temperature and we thus reject their temperature estimate. A 
white dwarf with a temperature of 13000 + 1 000 K would have 
been clearly detected in the OM, even at a distance of 140 pc. 

The unusual cool temperature of the white dwarf fits well 
in the current knowledge of all LARPs. Table |5] summarizes 
the white dwarf temperatures as determined so far. Apart from 



Fig. 11. Spectral fit of the combined PN and MOS spectra with 
a two component MEKAL model (see Tab. |4] for spectral pa- 
rameters) 



SDSS0837 all white dwarfs are cooler than llOOOK, which 
otherwise seems to be a lo wer limit for the white dwarfs in nor- 
mal CVs (' Gansickel2000l) . Without any prior accretion the tem- 
perature of the white dw arf in WX LMi impU es a cooling age 
in the order of 10^ yea rs ( Chabrier et al ]l2000l) as a lower limit. 
Schreiber & Gansickd (l2003b deduced a typical duration of the 
post common envelope phase of about 2 Gyr. Thus, CVs on the 
initiation of accretion due to Roche-lobe overflow should con- 
tain very cool white dwarfs. The low temperature also restricts 
the accretion rate in the past to have been less than 10 " Mq 
yr ', at least for peri ods shorter than the thermal timescale of 
the heated envelope (ITownsle y & Bild stenll2003h . Both limits 
support the assumption that WX LMi and the LARPs in gen- 
eral are pre-C Vs before the onset of Roche-lobe overflow. 

The UV light curves in Fig.|8]show a clear variability corre- 
lated with the visibility of the accretion spots, a clear indication 
for heated pole caps. Surprisingly the secondary spot seems to 
be hotter than the primary pole, which is inconsistent with the 
picture obtained from the optical spectra and the X-ray light 
curve. The reason for this is rather unclear. We can exclude a 
cyclotron contamination because the OM U-filter is centered 
between harmonics 4 and 5 of the second pole. 

The hardness ratio for the OM filters shows no significant 
variability which is somewhat confusing since our interpreta- 
tion of photometric variability of the UV requires also a hard- 
ness ratio variabihty. However, we cannot completely exclude 
that those expected count rate ratio variations are below our 
sensitivity threshold and cannot test this due to the lack of suit- 
able model atmospheres. 



3.4. The X-ray light curves and spectra 

The EPIC-PN light curve shows orbital phase-dependent mod- 
ulation with a pronounced minimum at the time, when both 
accretion spots are hidden from the observers view. While the 
cyclotron flux drops to zero around </> ~ 0.6, the X-ray flux re- 
mains at the constant minimum flux. We ascribe the residual 
flux around that phase to the corona of the secondary star. 
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Table 4. Spectral fits for the combined PN and MOS spectra 



Model 




NHP" 


kT (MEKAL) 


kT (BB) 


n'' 


Flux'' 








[keV] 


[eV] 


[cm-2] 


[ergs cm"- s"'] 


wabs (MEKAL) 


3.44 


1.3e-15 


0.62±0.02 




5.6e-6 


4.2E-14 


wabs (BB+MEKAL) 


1.42 


2.6e-2 


0.66±0.02 


117 + 11 


1.0'' 


4.6E-14 


wabs (MEKAL+MEKAL) 


0.98 


0.51 


0.26±0.02 




1.2e-12 


5.4E-14 








0.82+0.05 









" null hypothesis probability, * in units of le20, '" integrated flux 0.1-5 keV, '' frozen 



The average count rate of WX LMi is too small to make a 
phase-dependent spectral analysis feasible. Fig. [TTI shows the 
orbital mean spectra obtained with EPIC-PN and EPIC-MOS. 
Essentially all the X-ray flux emerges below 2 keV. We fit- 
ted the combined PN and MOS spectra with XSPEC. Our re- 
sults are summarized in Table |4] A first approach with a single 
MEKAL model for a coronal and/or accretion plasma, mold- 
ing emission from the secondary and from the accretion region 
with one temperature led to no satisfying fit. Motivated by the 
presence of a soft blackbody-like and a hard thermal compo- 
nent in high accretion rate polars, our next attempt included 
a black-body component. The fit was improved, but the tem- 
perature of 117 + 11 eV for the black body is higher than in 
any other polar observed so far and thus seems implausible. A 
satisfactory fit was achieved using a two component MEKAL 
model. The temperatures of 0.26 +0.02 keV and 0.82+0.05 keV 
are well within the regime that one could expect for a coronal 
plasma. The temperature thus determined would be very low 
for the accretion plasma of a normal-accreting polar It is how- 
ever in the same range as the temperatures derived from the 
analysis of the cyclotron emitting plasma. 

To make a possible distinction between the two components 
(accretion plasma and corona) we extracted spectra according 
to the time intervals when only one of the two spots is directed 
towards the observer The spectral fits of those two separated 
spectra yielded no significant differences between each other 
and compared to the overall mean spectrum. Also, we could not 
detect any significant variation in the X-ray hardness ratio as a 
function of orbital phase. Hence, we cannot discern between 
the coronal and the accretion plasma on the basis of possible 
spectral differences with the given number of photons. 

Figure [12] shows the spectral energy distribution from the 
near infrared to the X-ray regime for the individual spectral 
components. The graph shows AF^, i.e. the values plotted are 
representing the contributions to the total energy output from 
the system. Red squares indicate the contribution of the sec- 
ondary, which dominates the system absolutely. Black rhombs 
indicate the OM data and, together with the dashed line, indi- 
cate the contribution from the white dwarf, which can purely be 
observed in the ultraviolet. The spectrum shown in blue is the 
summed cyclotron spectrum from both poles. Arrows indicate 
the position of the cyclotron fundamental which are located 
in the //-band and between the H and K bands, respectively, 
i.e. the /T-band is likely uncontaminated by cyclotron radiation. 
The X-ray spectrum, shown as black solid line, is a mix of ther- 
mal plasma radiation from the accretion poles and the coronal 
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Fig. 12. Spectral energy distribution for WX LMi combining 
photometry in JHK (2MASS, red squares) representing the sec- 
ondary star, the OM onboard XMM-Newton (black rhombs) 
and the 8000 K nonmagnetic white dwarf model representing 
the white dwarf (0.6 Mq, d=100pc), the EPlC-PN spectrum as 
a mix of radiation from the accretion and the coronal plasma, 
and the sum of the two cyclotron components (blue line). The 
arrows indicate the wavelengths of the cyclotron fundamentals 
of the two poles. 

plasma. Accretion-induced radiation is clearly dominated by 
the cyclotron component. 

The unabsorbed flux (0.1 - 5 keV) for the two component 
MEKAL fit is 5.4 X 10"'''ergcm"^ s"'. For a distance of 100 pc 
this yields an X-ray luminosity of 6.4 x 10^^ erg s ' . Based on 
the variability of the X-ray light curves we assume that the sec- 
ondary and the accretion regions contribute roughly equal X- 
ray flux and luminosity to the observed X-ray spectrum, i.e. we 
get Lm„x ~ (Lix + i2,x) ~ 3.2 x lO^^ergs"'. 



3.4.1 . The secondary in X rays 

Our estimate of the X-ray luminosity of the secondary star 
makes a comparison with the X-ray activity of normal late-type 
dwar fs possible. The latter is correlated with the stellar rota- 



tion (Pall avicini et al.lll981l : iPizzolato et alj|2003h . According 



to Leggett et al.l ( 119961) the bolometric luminosity of our as- 
sumed M4.5 secondary is Lboi - 3.9 x 10^' erg s"'. The X-ray 
flux attributed to the secondary, Lm2,x ~ 3.2 x 10^^ ergs"', re- 
sults in Lx/L\,oi ^ 0.8 x 10"^. This result is only weakly depen- 
dent on the assumed spectral type, i.e. for an M3 secondary we 
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get 0.3 X IQ- 



Pizzolato et alJ (l2003h studied the relation between coronal X- 



ray emission and the rotation period. They confirmed the exis- 
tence of two regimes, one in which the rotation period is related 
to the X-ray emission and one in which a constant saturated X- 
ray to bolometric luminosity is attained (their diagrams 5 and 
6). This claim for the latest spectral types rests on rather small 
number of stars. The X-ray luminosity determined for the sec- 
ondary of WX LMi falls well on the constant branch of their 
diagrams, suggesting saturation of the activity. 

The same comparison can be made for the other low accre- 
tion rate polars as well. Table |5] compiles the parameters from 
the literature. Among the other systems only SDSS1553 has a 
clear X-ray detection and the source was found clearly above 
the Pizzolato et al. relation. However, the luminosity variations 



spectrophotometry of polars thus might further constrain the 
activity-rotation relation in a rather unexplored regime. This 
active secondary could be an explanation for the missing mag- 
netic white dwarf + main-sequence binaries in current catalogs 
jLiebert et al.ll2005h . If the magnetic white dwarf captures the 
wind from the active secondary, the emission of cyclotron ra- 
diation in the optical clearly vitiates the color selection criteria 
and/or the spectral composition. It requires a strongly magnetic 
white dwarf {B > 50 MG) to efficiently capture the wind, but 
the LARPs could mitigate the lack of those binaries at least for 
the case of highly magnetic white dwarfs. 

At the given orbital period, PqA - 2.78 h, one would ac- 
cording to Beuermann (2000) expect a secondary of spec- 
tral type M3 rather than M4.5. In this respect the secondary 
seems to behave as those in normal accreting pola rs which 



for the late spectral types in the compilation of Legget t et al.l also are observed typically 'too cool to comfort' ( Fri end et al 



19961) are large. Taking into account the uncertainty for the [l990l) . However, the standard explanation assuming the secon 



spectral type of the secondary, SDSS1553 could satisfy the 
Pizzolato et al. relation as well. Otherwise this could indicate 
additional X-ray emission from the accretion plasma, although 
the PN light curve shows n o variability, contrary to the optical 
Ught curve dSzkodv et all 12004). For SDSS1324 only^n up- 
perhmit was obtained from XMM observation dSzkodv et al.l 
2004 and for HS0922 there is an upper limit from a ROSAT 



PSPC observation (Reimers & Hagen 2000). SDSS2048 could 
be identified with a RASS s ource, but there are ot her optical ob- 
jects within the error circle CSchmidt et al ]l2005h . which makes 
the identification as an X-ray source ambiguous. If we take the 
upper limit as face value, it falls on the relation by Pizzolato. 
Hence, one (WX LMi) of the two objects with a proper X-ray 
detection seems to behave like a normal main-sequence star 
as far as X-ray activity is concerned and lends support to the 
picture of saturated activity, the other leads to no definite con- 
clusion, concerning the uncertainty in Lboi- 

4. Main results and discussion 

We have presented a photometric and spectroscopic analysis at 
optical, ultraviolet and X-ray wavelengths of WX LMi, a sys- 
tem termed in the past as a low-accretion rate polar (LARP). 
This designation was chosen due to its similarity to the polars, 
normal accreting magnetic CVs, as far as its stellar constituents 
and the magnetic field strength are concerned. 

Our main results can be summarized as such: WX LMi har- 
bours a highly active secondary star of spectral type dM4.5 
+0.5. The acitvity is obvious from photospheric Ha and chro- 
mospheric Ca ii lines, which vary in phase with the orbital mo- 
tion of the secondary, although not centered on the stellar disk. 
Activity is detected also at X-ray wavelengths at a level com- 
parable to other late-type stars at saturation. 

The phasing of the Ha-lines sheds some new light on 
old observations of polars in their low states. Faint residual 
Hcemission moving in phase with the s econdary (traced with 
the Nalines) was observed e.g. in U Z For ( Schwope et aljl997 ) 
or MR Ser jSchwope et alJll993h . although typically at a dif- 
ferent velocity than the Nalines. While Ha-emission was typ- 
ically interpreted in terms of irradiation, it may be that those 
low state polars revealed their active secondaries. Low-state 



daries being driven out of thermal equilibrium does not seem 
to work here, since the low temperature of the white dwarf, 
Tss < 8000 K, implies the absence of or a very low ac- 
cretion rate in the past. If WXLMi and its relatives would 
be just in an occasional low state of ac cretion the expe cted 
temperature should be above 11 000 K dOansickd I2OO0I) . as 
for the white dwarfs in normal CVs. The observed accretion 
rate is too small for a long-term equilibrium mass transfer 
dWebbink & Wickramasinghell2005l) . Also, the location of the 
second pole in WX LMi makes accretion via Roche-lobe over- 
flow rather unlikely. 



As suggested alr e ady bylSchwarz et alJ (1200 Ih for WX LMi 



and ISchmidt et alJ (l2005h for similar systems found in 
the SDSS the secondary is likely somewhat underfilling. 
Thus our analysis of WXLMi fully supports the view 
of LARPs as pre-Polars, post-common-envelope systems 
prior to Roche lobe contact, hk e a lready concluded b y 



Webbink & Wickramasin"gh3 (l2005h and lSchmidt et alJ (l2005h . 



We nevertheless detected ellipsoidal light modulations which 
in principle could be used to constrain the inclination and the 
filling factor, provided spectrophotometric data with higher ac- 
curacy could be obtained. 

We re-determine the distance to be = 100 + 20pc. 
The accretion luminosity, which is dominated by optical cy- 
clotron radiation with a minor contribution from X-ray ther- 
mal plasma emission, implies a low mass transfer rate of M ^ 
10"'^ M0 yr"'. Such low mass transfer rates are likely compat- 
ible with the wind mass loss rate from the a ctive secondary . 
Schmidt et alJ |2005), adopting a scenario by Li et al. ( 1995 ). 



proposed a 'magnetic siphon' channeling all material in the 
wind down to the pole caps with very little wind loss. The 
actual accretion geometry, determined by us from the orbital 
phases of the cyclotron features, doesn't seem to play an im- 
por tant role, once the magnetic field is sufficiently strong. 
Schwarz et al ] (I2OOI), modeling optical light curves and 



taking the beaming properties for the second spot into account, 
located the spots near the equator. We refine this model by a 
phase-dependent study of the cyclotron beaming for both spots. 
Primary and secondary accretion spot are located on the 'south- 
ern' hemisphere, i.e. away from the observer's hemisphere, in- 
dicating a field structure different from a centered dipole. The 
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Object 


Period 


Spectral type 


Comparison 


Distance 








Tv/D 


Reference"^ 










[pc] 




[erg/s] 


[10-3] 


[10^ K] 




WXLMi 


2.78 


M4.5 


G1268 


100 


3.86e31 


3.2e28'' 


0.8 


< 8 


4 


HS0922 


4.07 


M3.5 


G1494 


190 


1.97e32 


< 7.5e29 


<3.8 


< 10 


1,3 


SDSS0837 


3.18 


M5 


LHS377 


330 


3.43e30 






< 14 


3 


SDSS1324 


2.6 


M6 


GJllU 


450 


3.18e30 


< 1.4e29 


<46 


< 6 


1,2,3 


SDSS1553 


4.39 


M5 


LHS377 


130 


3.43e30 


1.4e29 


41 


~ 5.5 


1,2,3 


SDSS2048 


4.2 


M3 


G1388 


260 


8.93e31 


< 7.9e29 


<0.9 


< 7.5 


3 



" tak en one half of the un absorbed MEKAL+MEKAL flux in the range 0. 1 - 5 keV 
* see lLeggett et alj (Il996h 

' References: (l lWebbink & Wickramasinghd j2005l) , (2 lSzkodv et al.1 j2004h . (3 lSchmidt et alj j2005l) and references therein, 
(4) this paper 



field strengths in the two spots are 61.4MG and 69.6 MG, re- 
spectively. 

At such high field strength and low accretion rates as pre- 
dominating in LARPs probably no accretion shock forms. As 
shown by our deconvolution of the spectral energy distribution 
for the infared to the X-ray regime, the plasma cools mainly 
via cyclotron radiation instead of bremsstrahlung as in high- 
accretion r ate polars lending full support to the bombardement 
solution bv lFischer & BeuermannI ( 1200 lb . 

WX LMi as the brightest of the new class of pre-polars is 
a suitable target for spectropolarimetric observations at a 10m 
class telescope in order to search for polarized signal from the 
white dwarf's photosphere and to perform detailed modeling of 
the accretion plasma on its poles. 
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